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SUMMARY 

The Lighthill theory of aerodynamic sound requires a knowledge 
of the spatial and temporal variation of the two-point, two-time 
turbulent velocity correlations. The feasibility of determining 
these correlations based on extending closure models for one-point, 
one-time turbulence correlations is demonstrated. The procedure is 
based on a spatial moment integral formulation of the governing equa- 
tions using approximate, parameterized trial functions for the two- 
point, two-time velocity correlations. Solution of the equations 
results in a set of anisotropic length scales and the separation- time - 
dependent decorrelation of the ensemble averaged turbulent velocities. 
The analysis was simplified using the asstimption of homogeneous sta- 
tionary turbulence and a constant shear, unidirectional mean flow. 

It is shown that the anisotropic behavior of measured turbulence 
correlations can be characterized by this technique. Using the 
Proudman formulation of the Lighthill integral and the assumption of 
normal joint probability, measured sound power directivity can be re- 
produced for the compact acoustic limit by assigning a specific sep- 
aration-time behavior to a decorrelation function (which becomes the 
viscous dissipation in the limit of zero separation). 

It is concluded that the present approach is a viable technique 
for the prediction of turbulence generated aerodynamic noise. It is 
recommended that further effort should be concentrated on extending 
the theory to noncompact sound generation, developing a theory for 
the behavior of the decorrelation function, and to investigate more 
fully the effect of the anisotropic scales on sound generation. 
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1 . INTRODUCTION 


The theory of aerodynamic sound generation by turbulent flows , 
under certain simplifying assumptions, is based upon the avail- 
ability, either experimentally or theoretically, of two-point, two- 
time Eulerian velocity correlations. Theoretical analyses of the 
sound generation problem have generally adopted an isotropic form 
for these correlations. 

In addition to the form of the spatial and temporal variation 
of the velocity correlations, accurate predictions of sound genera- 
tion are dependent upon an ability to specify the one-point, one- 
time limit. The success of modeling efforts to close the turbulent 
rate equations for velocity correlations at second-order has pro- 
vided the means to predict these one-point, one-time correlations 
with greater confidence than previously. Bilanin and Hirsh^ have 
used such a model to predict the sound radiated from a turbulent 
swirling jet. Since only the one-point, one-time behavior of the 
velocity correlations could be predicted by their method, it was 
still necessary to adopt an assumed form for the characteristics of 
the velocity correlations in separation space and time. For this, 
the Ribner^ formulation was used. 

The feasibility of developing a more general representation of 
the two-point, two-time turbulent velocity correlations is explored 
in this study. If a Gaussian joint probability density function is 
assiamed, the turbulence statistics required to predict aerodynamic 
sound are the two-point, two-time turbulent velocity correlations 
Qi^(^,t,t,x) where t and t are absolute space and time and C 
ana t are the separation in space and time. Q.i(^,t,t,x) is de- 
fined by 


Q^j(x,t,t,x) = Uj^(x,t)uj(x + |,t + t) 

where the bar denotes ensemble time average. If the turbulence is 
assumed to be homogeneous and stationary, the velocity correlation 
functions are dependent only upon separation space and time. 

The equations governing Qij(x,|,t,T) contain higher order cor- 
relations in two-point, two-time triple velocity correlations, as 
well as correlations between turbulent pressure and velocity fluctu- 
ations. The closure problem is again the fundamental obstacle to 
the solution of the equations, as it is in the determination of 
single-point, single-time turbulent correlations. In recent years, 
however, methods of modeling the higher order terms in the one-point, 
one-time problem have been developed, resulting in significant im- 
provement in zero separation turbulent predictions compared to pre- 
vious mixing length models, as pointed out above. Successful 
application of these modeling techniques to the zero separation 
problem suggests that they may be extended to encompass two-point 
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correlations in space and possibly in time. 

Such an application presents numerous difficulties. The addi- 
tional variable in relative time- separation as well as relative 
spatial separation increase the difficulties in the solution to the 
governing equations many fold. Even this is no fundamental barrier 
to obtaining solutions, however. Higher order closure is again the 
basic problem, compounded by the additional spatial and temporal 
independent variables. 

In this study a new approach to the two-point, two-time turbu- 
lence problem has been attempted. Certain simplifying assiraptions _ 
have been made to reduce the complexity of the problem for this ini- 
tial study since testing the feasibility of the approach is the 
first priority. Therefore the turbulence is asstmed to be homogene- 
ous and stationary everywhere; conditions which later can be relaxed 
fairly simply in the absolute space variables under certain condi- 
tions. The governing differential equations are simplified by as- 
suming a constant unidirectional mean shear flow. The complexity of 
the equations is finally reduced to manageable proportions by employ- 
ing an integral technique. Approximate forms of the velocity cor- 
relations are chosen which are capable of reproducing the main 
features of experimentally observed correlations, each containing a 
set of separation"- time-dependent parameters which are essentially the 
anisotropic scaling factors. These approximate forms are substituted 
into the governing momentum equations and a family of moments in 
powers of the separation coordinates is taken over all space. The 
result is a set of simultaneous separation- time- dependent ordinary 
differential equations which may be solved for the time variation of 
the scale parameters. A set of subsidiary equations or constraints 
are employed which are derived from an integral form of the contin- 
uity equations. 

Once the scale parameters, as well as their time dependence, are 
found, enough information is available to determine the aerodynamic 
sound associated with this type of turbulent flow. Since the source 
term for the sound power is proportional* to the fourth derivative of 
separation time, a method must be devised to determine the fourth 
derivatives of the scale parameters. In this initial study we exam- 
ine the compact limit of sound generation and the theory will be 
applied at zero separation time. Thus, we avoid the necessity of ; 
integrating the governing differential equations and instead expand 
the parameters in powers of t , the separation time, and determine 
the coefficients of these expansions for small t . 

Although a number of simplifications have been introduced into 
the analysis during this initial study, they may be removed with ■ 
varying amounts of difficulty once the feasibility of the technique 
is demonstrated. Adoption of a working hypothesis that under special 
circumstances the turbulent behavior of a flow can be inferred from 
the analysis of simpler flows reduces this difficulty substantially. 
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For instance, if the characteristic size of the biggest eddies is 
smaller than the characteristic length of the flow over which the 
mean and turbulent variables may be considered homogeneous suggests 
that we may assume the flow to be homogeneous locally. Even if this 
condition is not precisely met, the homogeneous solution may serve 
as the first approximation to the actual solution. This point of 
view is adopted in this analysis. Thus, for instance, we shall as- 
sume that the velocity gradient is constant, and the solution for a 
varying gradient is obtained by a local application of the present 
analysis . 

This report is organized as follows: 

Section 2 is a synoptic of the approach taken. 

In Section 3 the equation defining the sound power inten- 
sity is introduced. 

In Section 4 the equations governing the two-point, two- 
time velocity correlations are derived. 

In Section 5 the models to be used for the pressure-velo- 
city triple- velocity , and dissipation terms are developed. 

Section 6 stmmarizes the one-point, one-time turbulence 
closure model used to determine the zero-separation 
velocity correlations which are the limit for the two- 
point correlations. 

In 7 the correlation function required to represent the 
approximate spatial and temporal variations of Q^j is 
selected. This is introduced into the governing equations 
in Section 8 and the appropriate equations are obtained. 

In 9 a qualitative comparison between theoretical results 
and experimental data is presented. 

In Section 10 the model is incorporated into the sound 
power integral and theoretical results calculated for an 
annular shear layer are compared with acoustic intensity 
data taken with an axis 3 mimetric jet in Section 11. 

In Section 12, conclusions and recommendations are offered. 


NOMENCLATURE 

A proportionality constant relating turbulent intensity 

to integral scale and shear, Eq. (6.7) 

A..,B.. directional factors in azimuthal (cl)) and latitudinal 
1JP9 (,|j) direction for sound power, Tables 5a and 5b 

A. . ,B. . directional factors averaged in (]) -direction in axisymmetric 
flow. Tables 6a and 6b 


4 



< p-l 


c(ipk> 
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gij(T) 

l(mk) 

ij 


1(4^) 


m 


N 

P 


P. . 
ij 

P(4Ji 4' .y) 
P(4^,y) 
(4j.y) 
q 

^13 

R. . 

IJ 

r 

s. . 

1-J 


low Reynolds nmber constant in isotropic microscale 
equation, Eq. (5.6) 

dissipation factor at high Re 3 molds number at zero 
separation; decorrelation parameter averaged over space 
for two-point, two-time viscous effects 

integral moments defined by Eq. (8.7) 

scale factor in direction 

ambient speed of sound 

J 

integral moments defined by Eq. (8.8) 
series expansion functions defined by Eq. (9.13) 
memory function for ij correlation 
integral moments defined by Eq. (8.6) 

total sound power intensity for axisymmetric flow per 
unit area in 4^ direction at observer distance x 

reference length 

v^/A 

mean pressure 
fluctuating pressure 

pressure-velocity correlations, defined by Eq. (4.19) 
sound power intensity per unit voltime 

(j) average of sound power intensity per unit of voltime 
turbulent intensity 

two-point, two-time turbulent velocity correlation 
defined by Eq. (8.4) 

vectorial distance from origin of separation coordinates 
axisymmetric jet radius to middle of annular shear layer 

triple-velocity correlations, defined by Eq. (4.20) 
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%(xj^ + x^) , absolute space coordinate 

T. . 

ij 

quadrupole strength density, Eq. (3.2) 

t 

absolute time 


mean velocity in k^ direction 


mean velocity in x^^ direction at 5]^ = 0 

■^0 

fluctuating velocity in k^ direction 


one-point, one-time turbulent velocity correlation 

e\ 

u, . 

IJ 

(u.u.) /q 

X j o'^ 

''f 

integration volume 

V 

c 

triple-velocity correlation modeling constant 

^k 

mean plus fluctuating velocity in k^ direction 


total velocity in direction of x 

w. . ,w. . 

ij ij pq 

defined by Eq. (10.7) 

X. . ,Y. . 

ij ij 

^k 

defined by Eq. (10. 13f) and Eq. (10.13g) 
absolute coordinate in k^ direction 

“ij ’^ij ’^ij 

correlation function separation time dependent 
scale parameters 

/\ /\ /S /S 

“ij ’^i.1 

°‘ij®ij^‘^l ’ ^ij'^ij^^a ’ ^ij°ij^3 ’ ^ij^ij ‘^l‘^2 

6. . 

ij 

Krone cker delta 

A 

s 

shear layer thickness 

r. . . 

ij 

dissipation of ij correlation 


defined by Eq. (4.5) 

K . . 

ijpq 

^ /N 

defined by Eq. (10.11) 

/s 

ij ijpq 

defined by Eq. (10.13c), Eq. (10.13d) and Eq. (10.13e) 
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isotropic form of integral length scale 

anisotropic integral length scale of ij 
correlation in direction 


turbulent microscale 


n(m^) 


viscosity 

kinematic viscosity 

separation coordinate; - x^) 

moment function defined by Eq. (8.5) 

density, ambient density 

decay length scale of ij correlation 

viscous compressive stress tensor 

separation time 

angular spherical coordinates 


Superscripts 


Subscripts 


vector quantity 
quantity at x' and t* 
quantity at x" and t" 
ensemble average 


reference quantity 

index of coefficients in Taylor series expansion 


value in k direction 
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2. SYNOPTIC OF THEORY 

The length and complexity of the theoretical turbulence and 
acoustic developments contained in the following pages make it de- 
sirable to summarize the approach more fully at this point. Lacking 
this, there is some danger that the objectives of the present study, 
and its important lines of approach, might be lost in the details of 
the analysis. 


2.1. Sound Intensity 

The Lighthill theory of aerodynamic sound^’*^ requires a know- 
ledge of the spatial and temporal variation of the two-point, two- 
time turbulent ensemble averaged velocity correlations Qij(^,t,t,x) . 
The sound intensity generated by the velocity fluctuation and radi- 
ated in the spherical coordinate direction ('p » <1>) » at distance x 
from the source, is 


I(ijj,(J)) 


1^; 2 5 
IbiT c 

o 


///».? 


X X 


dt 


( 2 . 1 ) 


V (x,t) and v' (x + t,t + t) are the velocity components at (x,t) 

and (x + 4,t + in the direction of the radiated power and con- 
tain the mean plus fluctuating turbulent velocities. The assximption 
of a normal joint probability distribution permits expressing 

2 2 

v^v^ in terms of , Q^jQpq ^iq^pj ’ other noncontri- 

buting terms. Obtaining Q.. is the objective of the turbulence 
analysis. ^ 


2.2. Governing Equations for 

Developing the governing equations for yields a set of 

differential equations having higher or der correlations uiuju]^ and 
pressure velocity cross- correlations , u^p ‘ . Equations tor the 
third- order terms contain terms in fourth- order velocity correlations. 
This is the so-called closure problem where equations for any order 
include correlations one order higher. An infinite hierarchy of 
equations result, each level of which does not contain sufficient 
information to obtain a solution. 


This problem is avoided by modeling the third-order and pressure- 
velocity correlations in terms of second-order correlations. Exten- 
sion of an established second- order closure technique is used for 
this, yielding a consistent set of (solvable) equations for . 

The equations governing the two-point, two-time correlations, 

, which are derived in this study, are 
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3Q.. 3Q.. 3U. 3U! , 3Q.. 

2 + Qjk 3^ + ^ik 3l^ + + ^k^ 3^ 


where 


3Q. . 

+ - “k' - 7 


I sir 


lUiUjUk 


^ ■ IF 


3f: ^jP + 3iT 

1 J 


— u.p' 

p 1^ 3^j 


3q 


1 3^Q.. 

+ -^v — + 2v 

3s^ 




( 2 . 2 ) 


^k ' ’^k - ’'k 


®k 2^^\ ■*■ ’‘k' 


T = t* - t 


(2.3a) 

(2.3b) 

(2.3c) 


Assuming that the turbulence is homogeneous and steady eliminates 
derivatives in s]^ and t and assunption of a constant, unidirec- 
tional mean velocity gradient in the X 2 direction 



(2.4) 


affords considerable simplification. The following turbulent models 
are then introduced. For the triple-velocity correlations 


3 

3 5 k 



u. 1 




(2.5) 


The pressure velocity correlations are modeled as Rotta tendency-^ 
towards- isotropy terms 
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TTTr- U.p' - 
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= £ 
A 


[q^ - 


( 2 . 6 ) 


The dissipation model (high Reynolds ntimber limit) is 


2w2q.^ = -2b a Qy 


(2.7) 


Here A is a global macroscale set by the local flow and q = u.u. 
The equations then take the form ^ ^ 


3Qf. / 9Ui \ ^ / 




jdX2 


" a . - ^6..Q 

m 


ij 


+ V qAV^Q. . - Q- • 
c^ A ^xj 


( 2 . 8 ) 


b and are modeling constants. The form of the modeled terms 

are adapted from Donaldson's second-order closure technique.^ 

2.3. One-Point, One-Time Correlations 

The limit of Qij for C and t equal to zero are the one- 
point , one- time velocity correlations. The equations governing 
Qii ( 0 1 0) may be solved to determine the one-point, one- time values 
or u^u. in a constant shear flow. The results are 


^lUi/q^ = -j[l + Ab] 




^3^3 1 

= |[1 - 2b] 


2 1 


UiU2/q 3A 


where A = q/ (AdUj^/dx 2 ) . 


r 3^ 

- 1 - 2b l 

= [ 3b J 


(2.9a) 

(2.9b) 

(2.9c) 

(2.9d) 
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The one-point, one-time solution used in this study is for con- 
venience only. The one-point, one-time correlation can be obtained 
from any turbulence model, or from numerical solutions of shear flow 
turbulence. 


2.4. Two-Point, Two-Time Correlation Functions 

This study employs an integral approach to determine the behavior 
of Qij(?k'^) • analytical form is selected for Q^j which is 

capable of approximating the two-point, two-time behavior of velocity 
correlations seen in available experimental data. This form contains 
a set of separation- time- dependent free parameters which are deter- 
mined using the governing differential equations together with an 
integral approach. Here, we have chosen the following form to approx- 
imate 



u . . R . . g . . ( x) 



[1 - 


a . . 

xj '^1 


- U^T) - 




- Y • -?o 
xj^3 


- ^ij<^l 


- U^T) 


^2^ 




( 2 . 10 ) 


where u^Uj is the zero separation, one-point, one-time result. 

contains the parameters a^- , Bj[j , Yi-j , Pij , <^ij > cq , C 2 
ana c-^ which are separation time dependent." a^j , , yii and 

yj_j serve to provide zero crossings of the correlations, the ani- 
sotropic analog of the isotropic velocity correlation function g(r) 
zero crossing. The asymmetry in the ^2 direction due to shear is 
accounted for by the presence of Uij . ^^,02 and C 3 are an- 
isotropic scaling functions to provide variations in spread in the 
three coordinate directions. Sij^'^^ is a memory function which 
decorrelates Qij at ?k “ 0 in a coordinate system following the. 
flow at velocity . In order to evaluate the anisotropic be- 

havior of the turbulent scales , we define 


/ ^ij^^k ■ '^Ik^l t) , k = 1,2,3 


3-J 


’ij 


7 


k ^Ik^l 


o (2.11) 

which becomes a function of the correlation function parameters. 
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2.5. Continuity 


The continuity equations governing 




0 


( 2 . 12 ) 


These equations serve to add constraints on the parameters and pro- 
vide continuity in an integral sense as will be demonstrated shortly. 

2.6. Integral Moments 


The correlation functions are substituted into the governing 
equations and a set of integral moments of the equations are taken 
over all space. The moment function 


l(mk) 

ij 


OO OO 00 

-III’' 

— 00 — 00 — 00 


(mk)Qij 




(2.13) 


is defined where it (mk) is a function of the r coordinates 

k 


H (m^ yXTL^ »m^ ) 


^mi^m2^m3 


(2.14) 


A sufficient mmiber of moments are taken to provide enough equations, 
together with the continuity constraints , to allow solution of the 
set for all the unknown free parameters. The governing equations 
become total differential equations in t in the separation time 
dependent parameters. Since experimental data indicates that 
3 QjL 4 / 9 T = 0 at zero separation time, the initial conditions are 
found by an algebraic solution of the equation set for the parameters, 
with zero rate of change at t = 0 . The characteristic length of 
the flow is specified by independent selection of one scale of the 


set • 


Continuity is satisfied in an integral sense by integrating 
the continuity equation over half the separation space domain as 
follows 



pqr = 123,231,312 


(2.15) 
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2.7. Decorrelation Function b(t) 


Decorrelation of Q^j with separation time is provided for in 
this study by a separation- time dependence of the turbulent dissipa- 
tion function b(x) . At zero separation in space and time, b is 
a modeling constant. The one-point, one-time value of b = 1/8 is 
probably no longer adequate when considering a space integrated dis- 
sipation of Qij . Instead a time varying form b(t) is hypothe- 
sized. b(0) and the time dependence is selected to provide 
agreement between the theoretically calculated sound generation and 
measured sound intensity data. 

2.8. Solution of Time Dependent Equations 

The resulting equation set may be solved ntmierically . 

However, for this first effort the compact acoustic limit 
is assumed. This permits expansion of the parameters in power series 
in T about t = 0 . The zero-order solution is the initial con- 
dition now known. Higher-order coefficients are found by substitu- 
tion of the expansions in the differential equations and ordering of 
terms. The result is a set of algebraic equations for coefficients 
of order n + 1 in terms of coefficients of order n . The decor- 
relation function b(x) is specified by selection of b(0) and the 
higher order coefficients to provide agreement between measured and 
calculated sound intensity distributions. 

2.9. Calculation of Sound Intensity 


Having the behavior of Qi* ■? for the selected shear flow con- 

2 2 

sidered, the results can be used to evaluate v^v^ (^lu- , 

dUn 


XI . X2 , X3 


dx2 

with respect to t . Then t 
tion taken over separation space. Since 


, (J)) and the fourth-derivative is then taken 
Then t is set equal to zero and the integra- 




is an explicit func- 


tion of this is done analytically. We thus have the compact 

limit of the locally radiated sound intensity, which is a function 
of the lo cal gradient dU]^/dx 2 and the one-point turbulence corre- 
lations U£U4 . A final ntmerical integration over a specified mean 
velocity profile and turbulence distribution provides the sound inten- 
sity radiated in the direction. 

This has been a very brief overview of the theoretical approach. 
The analysis is presented in considerably more depth in Sections 3 
to 11. 
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3 . 0 ACOUSTIC SOURCE TERM MODEL 

The sound pressure radiated to a point x in the far field in 
a localized unsteady or turbulent flow was shown by Lighthill to be 
given by 


p(x,t) = ■ 


^ f\ 

2 

■3^T 

O CO 

. 3t 


il 


dy 


T.. is a quadrupole strength density 


T. . = pv.v. + T . . + (P - 
ij 1 J ij 


c2p) 


6 . . 


(3.1) 


(3.2) 


where Pv^vj is the unsteady momentum flux, vi the velocity, t.. 
the viscous*^ compressive stress tensor, P the local pressure, 
the ambient speed of sound, p the density, and the Kronecker 

delta. Equation (3.2) is normally dominated by the unsteady momenttcn 
flux pv^v^ . The symbol [ ] in Eq. (3.1) denotes retarded time. 

The indices i, j and k are equal to 1 , 2 or 3 and repeated indices 
are sxjmmed. Figure 1 presents a schematic diagram of the coordinates 
employed here. 

The sound power generated by the velocity fluctuations and radi- 
ated in the directions (i^,*!)) in spherical coordinates is the ensemble 
averaged pressure fluctuations at the point of observation divided 
by Po^o • Defining this by I(ip,<fi) 


I(ifj,(t)) = 


XiX.Xj^Xi 

77 2 T~B 

l^TT p^C^X Ob "00 


// 


3^(pV^vM 9^(p’'vj;v’{) 

2 

3t^ 3t 


t' - t" is the difference in time of travel to 
y" . Ribner® expressed Eq. (3.3) as a function 
and the separation in space and time using 

y = |(y' + y") 

t = y' - y" 


T = t' - t' 


dy’dy” 

(3.3) 

at y',f 

and 

5le average. 


from y' 

and_^ 

the midpoint y 


(3.4a) 


(3.4b) 


(3.4c) 
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Figure 1. Coordinates of acoustic analysis 
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and the assumption that the observer distance x is large compared 
to the flow dimensions, i.e., 


c^T s ^*x/x 


(3.5) 


giving 


x<rhere 


■I 






(3.6) 


P (<f» , <f) »y) = 

leTT^c^x*^ j 3 t^ 1- J K i 


O 00 


Proudman’ expressed Eq. (3.7) in the very convenient form 


P(i|j,<{> »y) 


16 fr‘ CoX 


,5. 

-o" 


■f— 

■J 


2 

V V 
X X 


dt 


(3.8) 


where Vx^ and are the components of the velocity fluctuation at 

y' and y" in the direction of x . 

The correlation v^v^ can be expressed in terms of the quad- 

rupole correlations v^j^VjV^V]^ and the coordinates x,ij^,<|) . Assump- 
tion of a normal joint probability distribution for u^ and u^ 
permits expressing the fourth-order correlation in the terms of 

Q.. and 6, , . This will be done in Section 10. 

^kl 

4. DEVELOPMENT OF TWO-POINT, TWO-TIME, 

TURBULENT CORRELATION EQUATIONS 

Here we develop the equations governing the ensemble averaged 
two-point, two-time correlations of turbulent velocity components 
that are required to predict aerodynamic sound. As we shall see the 
equations for the two-point, space- time velocity correlations con- 
tain triple-velocity and pressure- velocity correlations which cannot 
be solved for exactly without recourse to an infinite hierarchy of 
equations, each containing successively higher-order correlatxons . 
This is the so-called closure problem. It will, therefore, be neces- 
sary to model these terms in a way which will meet certain criteria. 
Now, in the limit of zero separation in space and time the exact 
equations governing the correlations reduce to the equations des- 
cribing the Reynolds stresses at a point. We shall model the two- 
point, two-time equations such that at zero separation and time they 
reducG to an established theory of higher-order closure for the 
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Reynolds stresses. Such closure theories have been developed by 
Donaldson®, Hanjalic and Launder®, Wolfshtein, Naot and Lin® and 
Lumley and Khaj eh-Nouri^ ° , among others. Closure theories generally 
specify numerical values of coefficients in their modeled terms 
chosen to optimize agreement between theory and experimental data. 

It is considered undesirable to require changes in these coefficients 
for various flow conditions. The coefficients in our modeled terms, 
which as we shall see are the viscous, pressure-velocity and triple- 
velocity correlations, will be assigned those values selected by^ 
Donaldson and his associates and employed in their technique of in- 
variant modeling. ® ’ (Henceforth this technique will be referred 
to as (I) for convenience. ) 


4.1. One-Point, One-Time Equations 

The equations governing the conservation of momentum and mass 
in an incompressible uniform density fluid with constant viscosity 
are 


8U. 3U. 

i + U 

at ^ 3x^ 


1 ^ 
P 3x^ 


+ V 


9x? 


3U. 
1 

3x. 

X 


= 0 


(4.1) 


(4.2) 


The velocity Ui is now written as the sum of mean and fluctuating 
components, + u-j_ , introduced into the governing equations, and 
an ensemble average is taken. The results are 


3U. 3U. 

— - + U — - = 
3t ^ 3Xj 


3u.u. 

3Xj 


^^i 3x? 


(4.3)- 


9U. 

X 


(4.4) 


These are the equations governing conservation of the mean mass and 
momentum at one-point and one-time. The equations for the Reynolds 
stress are developed from Eq. (4.1) and Eq. (4.2) by introducing the 
mean and fluctuating components for the velocities u^ and Uj and 
before ensemble averaging, multiplying the respective equations by 
u. and u. . Addition of the equations, subtracting out the mean 
u^ing Eq. U-3) and Eq. (4.4) and averaging l ead to the equations 
for the ensemble averaged Reynolds stresses Uj^Uj at a point . 
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U . U, -:r— ^ 

J k 3x^ 


3u .u. 
X 1 

3t 




3u.u. 

3-JL 

3xi. 


3U. 




U. 

1 

P 3x. 


(4.5) 


J. 9£_ 

P 9x^ 


+ V 


2 

3 u . u . 
^ J- 

3x? 


-2v 


3u. 
: 

3x, 



These equations are the counterparts of the equations governing 
the evolution of the two-point, two-time velocity correlations, which 
will now be derived in more detail. Initially we do not assume homo- 
geneous turbulence, but for simplicity this assxraiption will be made 
later, as well as that of a constant gradient mean flow along one of 
the coordinate axes. 


4.2. Two-Point, Two-Time Equations 

Now define the mean and instantaneous turbulent fluctuating veloc 
iti0s at point xj^ and time t as ^i, s^ud ^i, » sud their counter- 
parts at point x^ and time t' as Uj and Uj , introduce them in 
to Eq. (4.1) and subtract out the mean equations. This leaves 


u. 3u. 3U. 3u. 

t“ 3x^ ^ 3xJ^ ^k 3x^ 


3u_. 
3 




i l£_ 

P 9x^ 


+ V 


3^u, 

3 Xu 


(4.6) 


9u'. 


aui 


3UI 


9u'. 


3t' ^k 3x,'_ 3x^ ^ 3 


xV 




= .- i 14 + V 

p 3Xj 


3^ul 


9x. 


,2 


(4.7) 


Now multiply Eq. (4.6) by ul and Eq. (4.7) by u^ . Sxnce primed 
dependent variables are not functions of unprimed independent varx- 
ables, and vice versa, they pass through the partial derivatives. 
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Adding the two equations and taking the ensemble average gives 


Bu.ul 3u.u! 


9t 


^ + U. 


9u.u! 9u.ul 9Uj 

_1_JL 4- n' — 4- ii'i 


9t' 


k 9Xr 


\ 9x/^ 9xJ 


9Ul 9u.u,ul 9u..ulu,' 

^i^ Tx?" W, 


iVi ^ 

“ 9x,' 


(4.8) 


1 9 


P - 


u:p - 


— U.p' + V 

p 9Xj 


9 '^u . u I 

3. 1 

9x? 




9x. 


,2 


New independent variables are now introduced to differentiate 
between the effects of absolute position and time and separation 
distance and time. Define 


‘’k 

®k 

T 


^k ■ ^k 


2^^k ^k^ 


= t' 


(4.9a) 
(4.9b) 
(4 . 9c) 


The derivatives become 


9 

9Xk 

3 + 

3Xk 3Sj^ 3X^ 

9 _ 1 

9?k ■ 2 

9 9 

9Sk 9^j^ 

(4.10a) 

9 

8 x^ 

*®k 3 ®^k 

3X^ 3S^ 3X^ 

9 _ 1 

3Ck 2 

9,9 
9Sk 9Ci, 

(4.10b) 

_9i 

- 1 3^ y .3^. _ 

32 


(4.10c) 


“42^ 2 

9s^ 9?k 

38^3 Ck 



_ 1 9^ . 9_^ ^ 

32 


(4.10d) 

3X^2 

“4 2 2 

9s2 95^ 

9Sk9?k 
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3 t 3 ^ 3t 3 


3 ^ ^ 


3t’ 3t' 3 t 3t' 3t 3 t 3t 


(A.lOe) 


Taking the ensemble average we arrive at 


3u.u'. 3u.ul 


3U. 3U! 

j“k ?ir + “i“k 15+ 


3u.ul 9u.ul 

+ 1 t“k + Uk' -st:^ + i«k - 


- 1 sir '“i“j“k + “A"!' - si: ‘'^i“j“k - “i"k“i' 


1 

2p 


3J-- ’^jp + rir "iP'l ■ p 


35^ “i 


u,p' - — 


1 .. 3- 


+ -O V 


uTuT + 2v 


3s 


2 LA • . 

1 J 
k 


u.u'. 


3C 


2 1 J 


(A. 11) 


A. 3. Simplifying Assumptions 

We now invoke the assumption of homogeneous and stationary tur- 
bulence, and all derivatives with respect to the absolute coordinates 
a nd absolute time t vanish. Defining 

u^(x^ ,X 2 »X 2 , t)uj ^1 ’ ^2 ^ ^2 ’ ^3 ^3 ’ ^ ^ 

we obtain 

3Q.. 3U. 3U! 3Q.. 

rF- + + Q,-v ^k^ 


3t 


^jk 33^; ^ ^ik-^ ^"k '^k^ 


= - sir t^i^K ■ 


1 

P 


917 “iP 


_3 

3^ 


rip] 


+ 2v 2 ^ij 


(A. 12) 
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Figure 2 summarizes the coordinates systems in absolute and 
separation coordinates. 


One of our basic working hypotheses is that the behavior of more 
complex flows may be inferred from the analysis of simpler flows, 
specifically here that locally the turbulent correlations can be cal- 
culated from knowledge of the local velocity gradient. This gradient 
is taken to b e con stant over the volume of separation coordinates 
within which u^uj becomes completely uncorrelated. The convective 
velocity in Eq. (4.12) is assumed to be along the coordinate 

with a constant gradient in the direction. Then 



U 


k 


so Eq. (4.12) becomes 



dU^ 

^lk^2 dx2 


(4.13) 

(4.14) 


3Q. . 

JJ. + 


3T^ + <'^11^23 




dUj 

dx. 


+ K 


dU, 3Q.. 
1 3-1 

2 dx2 3^2^ 









+ 2v 




(4.15) 


This equation describes Qij with respect to a coordinate system 
convected with the local velocity within the shearing flow. If 
we wish to write the equations for a flow moving at a velocity Ui^ , 
with respect to the coordinate system, then we have 


3Q . . dU, / dU, \ 3Q . . 


H, 






(4.16) 
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Now, for homogeneous turbulence in the coordinate system adopted 
and since the turbulence is stationary, the correlations will not be 
changed if we replace by -K, and t by -t . From the con- 

dition of invariance under translation, for any homogeneous flow 
field* 


uTu^uT (5^ ^ 

u^u^Uj ~ ^i^k^j 

pUj (^2^>^2’^3*^^ ~ P 

From the condition of invariance under reflection for this co- 
ordinate system 

u . u j 5 2 ’ ^3 ’ ^ ^ ~ ^i^j ( “ C-j^ » ~ ? 2 ’ ~ ^3 * ^ 

Uji^u^Uj (?^i^ 2’^3’'^^ ~ ~ ^i^^k^j ( ~ “ C 2 > ~ ^3 ’ ~ ^ 


u^Uj^Uj ~ ^i^k^j 

pu! (5^,52>^3>^^ ~ ~ (“^2>“^2’~^3*”"^^ 

The triple-velocity correlations are S 3 mmietric with respect to the 
indices referring to the same point so that 



^k^i^j 





~ ^ ^ 2 , 2 
*The following symmetry arguments can be found for t = 0 in Hinze , 

p. 332. 
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The term in Eq. (4.17) containing the triple-velocity correla- 
tion term gives 

3 [u.ulu7 - u .u, ul ] g =0 

^ T = 0 

From continuity, contraction of the pres sure- velocity correla- 
tion term gives 


“ h 7 ^iP ■ ° 

Contraction of Eq. (4.16) then yields 


9Q dU, 

-57- + 2Qi2 


/ dUi V 3Q . . 

(\ * ^2 


9^Q. . 


1C d^y. 


(4.17) 


2 2 

Now at zero separation 9 written 


3 

9C 


2 ^ii 
k 


3x^3x^ ^ 


3u. 3ui 

53^ sSJ 


SO that for = 0 ,t = 0 , Eq. (4.17) reduces to 


3u.u. dU, 3u.u. 

3^ ^^1^2 ^1^ 3q“ " 


3u. 3u. 


= - 2v 


3xj^ 3xj^ 


Since 5 


dU. 3u, 3u. 

1 = _ i i 

1“2 dx2 9x^ 3x^ 


Ui u 


(4.18) 


at = 0 , T = 0 . 

To facilitate writing equations let 


*■11 - [t|t “iP' - sir 


(4.19) 
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Written out in full Eqs. (4.16) now appear as 
8Q, , / dUi \ 3Qi 1 dU 


(4.20) 


^^22 / 


^ +fu 

Vlo 


+ 


^Qi2 / 

^ <“lo 


+ 



> ''^11 _ 

dX2^ 


dUjN 

1 *'522 - 

dX2y 

1 

dUi^ 

1 '?33 _ 

dx^; 

’ 'h 

dUiN 

\ ^^12 

dX2, 

/ 9^1 


!i2 35jr- ^11 - ifii + 


- ^22 - 7^22+ 


- s 


33 


dUn 


= - Q 


S,o 


22 dx2 P 


ip 33 + 2 vV^Q 33 ( 4 . 21 c) 


- ip32 + 2 vV^Q^ 2 (4.21d) 


Although equations may be written for the Q]^3 and Q23 correla- 
tions (which are equal to zero at zero separation), they may be found 
from continuity, given the remaining correlations. We shall neglect 
Q13 and Q23 in this initial study for simplicity, concentrating 
only on the correlations which are explicitly required to compute the 
energy components Qn , Q22 Q33 • 

Although it is not possible to define precisely the character- 
istic physical behavior of each term in the governing equations for 
the two-point, two-time velocity correlations, they may be identified 
by analogy with the one-point, one-time equations as follows. Write 
the equations as 




^ " (“lo 


2 dx2; 351 


(«ilQ2j + ^ji'^12) dxj 


I 






> 


2 


“tl 


'm' 


+ 2 


V 


2 Qij 

IV 


(4.22) 
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I is a production- like term, whe re in teract ion w it-h the mean 
flow generates a net increase in the u^ui and uiU 2 correlations. 

II, the pressure-velocity correlations, disappear upon contraction. 

This term transfers correlations between components. In one-point, 
one-time theory it is the "tendency -towards -isotropy” term. Since 
production of energy -containing eddies occurs in the i=j=l equa- 
tion, this term transfers energy to the 22 and 33 components which 
in turn may exchange ^energy among themselves and the 12 correlations. 
Deissler^^ and Fox^^ have studied this transfer for low Reynolds 
number turbulence. A summary of their results appears in Hinze^^ . 

Ill is the triple-velocity correlation and is a diffusion- like 
transfer of Q^. by turbulence gradients. Finally, IV is a dif- 
fusion-like term which serves to decorrelate through the action 

of viscous stresses. At zero separation it redunes to the dissipation, 
9u, 9u. 

-2v-r— ^ -r— ^ , and the triple- velocity correlation term is zero. Thus, 

° ^k ° ^k 

the contracted form of the equation, for which the pressure-velocity 
terms sum to zero, shows that the dissipation must be balanced by the 
production for stationary turbulence. 

5. MODELING OF PRESSURE-VELOCITY, TRIPLE-VELOCITY 
AND DISSIPATION CORRELATIONS AT SECOND-ORDER 

In this initial study, intended to test the overall validity of. 
our approach, it is desirable to provide closure of the governing 
equations using the simplest possible model consistent with physical 
principles. The philosophy here is to extend the form of one-point, 
one-time modeling to the two-point, two-time correlations such that 
the latter reduce to the former for = 0 , t = 0 . 

5.1. Modeling of Triple-Velocity Correlations 

The triple-velocity correlations at zero separation in space 
and time have usually been modeled as gradient diffusion terms in 
previous closure techniques. The form adopted here and which has 
been used in (I) is 


9 




(5.1) 


where = Uj_U£ and A is a turbulent length scale. In anisotropic 
turbulence we should expect that A is dependent upon direction. 
Indeed, this would be the natural extent ion of Eq. (5.1) to a more 
sophisticated model. Later, we shall define a set of directionally 
dependent turbulent scales. In the present turbulent diffusion model 
A will be assumed independent of direction, although we might expect 
A to be dominated by the assvmtied unidirectional mean velocity gra- 
dient in the X 2 direction. 
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Equation (5.1) satisfies the ten sor S 3 raimetry of the triple- 
velocity derivatives 3 /9xk(uiujuk) but not the symmetry of u^u^uk . 
The model coefficient v^ was assigned the value of 0.3 in 
(I). Th is, mcp del is extended to a two-point, two-time modeling by 
writing njuj in place of the one-point, one-time correlation. Due 
to homogeneity the turbulent diffusion coefficient v^,qA is assumed 
uniform in space and so 


Sij = - v^qAV^Q^j (5.2) 

where 



5.2. Modeling of Pressure-Velocity Correlations 

As noted previously, contraction of the governing equations re- 
sults in elimination of the pressure~velocity correlations since 
Pii = 0 in homogeneous turbulence. For one-point, one-time correla- 
tions these terms redistribute energy between the velocity components, 
tending to decrease anisotropy and decrease the turbulent shear 
stress. According to Deissler and Fox, P^. may actually promote 
anisotropy in two-point correlations under^some circiamstances . 

From the governing equations for i = j we see that the pro- 
duction term appears only in the equation. Py must therefore 

increase Q22 Q33 the expense of Qll . Rotta^ ® modeled 

the pressure-velocity correlations as a tendency- towards -isotropy 
term for one-point, one-time calculations as 


1 

P 


u 


9P 

i3Xj 


+ 





u.u. 
1 J 


- 


where C is an order one constant. Taking C = 1 
Eq. (5.3) is extended for the present analysis to 


(5.3) 


the form of 


which reduces to Rotta's form at zero separation. Note that Qij 
could be transferred unequally among components by writing Eq. (5.4) 
as 


-V . . 

P 


= SL 
A 






(5.4a) 
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where we must set =1 to make = 0 . For simplicity in 

this initial study we shall use Eq. (5.4). 


5.3. Modeling of the Dissipation Term 

2 0 

Term IV of Eq. (4.22), +2v° Qij/95k » is a diffusion- like trans- 
fer of correlations due to visc o us effects. At one-point, and one- 

C) 4 S'Xln 

time this terms becomes under contraction of the 

indices balances total production. For large Reynolds number the 
magnitude of this dissipation forces a high rate of decorrelation in 
near zero separation, and therefore the second derivatives of 
Qij are large in this region. If we express Qij near zero separ- 
ation as 



(u.u. ) 

1 J O 


^i^i 


+ 



(5.5) 


where X is th e turb ulent microscale, then the dissipation is pro- 
portional to -(u£Uj)o/x2 • As suggested by Rotta^ ^ ^ ® , the micro- 
scale is modeled as 


X 


2 



a + 


bqA 

V 


(5.6) 


Therefore, in the high Reynolds number limit the dissipation term 
will be 





(5.7) 


where (u^u^ ) has been generalized to Qy . The one-point, one- 
time limit of the contracted form of this model is -2bq3/A . 
Donaldson in (I) has used the isotropic form of this model, (i.e., 
2/36j.bq3/A) , with b = 0.125. For non-zero spatial separation we 
might'^ expect a smaller value of b since the viscous decorrelation 
occurs over a larger characteristic length than X . For increasing 
separation in time it is reasonable to suppose that decorrelation 
will increase. For this study we shall assume that b is some aver- 
age value over separation space and a function of t . A prelimin- 
ary value will be assigned at some point as well as coefficients in 
the separation- time expansion. These values will be chosen to pro- 
vide physically realizable behavior of the correlations and, to match 
experimental acoustic measurements. 
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5.4. Modeled Form of Equations 


Equations (4.22) now take the form 


/ dU, \ 8Q . . / 

(”io ' "(*ii 


2j + '^jl^i2y dx2 


iWj 


Written out in full these are 


+ Al + P y^ll - 20 ^ [20 -Q -Q 

TT (^^1^ ■*■ ^2 d^ya^i " ■ ^^12 dx2 3A [^^11 ^22 ^33 


h2+(^ +5 

3^ + 4 dx2;35l 


+ v^qAV^Q^^ - ^ 


■ 3A [^^22 ^ll'^33 


+ v^qAV^Q22 " ^ Q 22 


!i3 ^ 

,T M\^^2dx2hCi 


“ ■ 3A ^^33"^11'^22] 


+ v^qAV^Q 33 - ^ Q 33 


^12^/u 

)T ^2dx2)9?^ 


dU^ ^ 

^22 dx2 ” A ^12 


+ v^qAV^Q^2 ■ ^ ^12 
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6. ONE-POINT, ONE-TIME TURBULENT VELOCITY 

CORRELATIONS FOR STATIONARY, HOMOGENEOUS 
SIMPLE SHEAR FLOW 

In applying the present theory to calculate turbulent acoustics 
the values of q and A in Eqs. (5.9) would be obtained from full 
nijmerical solutions of the flow being studied. For convenience, and 
to provide a consistent set of one-point, one-time velocity correla- 
tions that will be required in order to evaluate the shear flow model 
we are developing here, l et us examine the limiting form of the 
governing equations for (uiUj ) at = 0 > = 0 . These correla- 

tions are the initial conditions for the two-point, tx<7o-time correla- 
tions. What follows is adapted from Donaldson's superequilibriun theory in (I). 

Since the flow is homogeneous and stationary, derivatives of 
(uj[uOQ with respect to absolute space and time are zero. At the 
limit of zero separation in space and time, Eqs. (4.5) define the 
one-point, one-time turbulent velocity correlations for this flow. 

The left hand sides of the equations vanish. The viscous term trans- 
forms to a dissipation term at zero separation, dependent upon the 
curvature in separation space,' of the two-point velocity correlations 
at zero separation. The equations at zero separation become 


dU, 


0 = - (6 


.,u„u. + 6.,u.u«) j— 
rl 2 j jl 1 2 dx. 


- ^ [u.u. - 
A 3- J 


n2l 

Q o • • Q J 
3 


- 2v 


3u^ 3Uj 


( 6 . 1 ) 


where the modeled form of the pressure-velocity correlations is used. 
The contracted form of this equation yields 


dUj^ 3u^ 3u^ 

^ 1^2 


( 6 . 2 ) 


or production exactly balances dissipation, where e is the total 
dissipation rate of the turbulent kinetic energy per unit mass, 
l/2q2 . If we define the dissipation terms in general as 



3u^ 3Uj 


(6.3) 
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then the equations may be written 


2e - ^ 
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■ ^^12 " 

0 
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, 4 d) 

contracted equations 

show, of course, that 


- 


e 

== e . • 
11 

Ell + 

^22 

33 


(6, 

. 5 ) 


Equations ( 6 . 4 a), ( 6 . 4 b) and ( 6 . 4 c) are, therefore, not independent 
since any one can be derived from the other two. 

Manipulation of Eqs. ( 6 . 4 ) yields the velocity correlations in 
terms of the dissipation 


UfUi 


1 , 2A (p 4. p 'k 

3 3 '■ 22 33^ 
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(6 . 6a) 
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(6.6b) 
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dUi 


dx. 


2A 


.3 ^22 
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(6.6c) 


(6.6d) 


We can solve for £^2 terms of the other dissipation components. 
First let us define’^ 
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(6.7) 


A = 


A 


"lu^ 

dx. 


which for a constant value of A is a definition of A . The ratio 
q/dUi/dx2) can be considered a characteristic length in a constant 
shear turbulent flow. The inverse of A then defines the size of a 
"typical” eddy in terms of this characteristic length. 


e ^2 then becomes 


12 


= ialfcl 

A A [ 4 


A^ E + E ) 

2 a. e t- ^3 


-i] 


( 6 . 8 ) 


In shear flow we would expect to be anisotropic. Past 

investigation of shear flow have usually assxmed isotropic dissipa- 
tion 


e . . = -^ 6 . • b ^ — 
ij 3 ij A 


(6.9) 


For simplicity we shall also use this assiraiption and Eqs. (6.6) give 


UiUi , 

= ^[1 + 4b] 


(6.10a) 


^2^2 


U3U3 


^ = f[l - 2b] 


^ l ’^2 


3A 


[1 - 2b] 


(6.10b,c) 

(6.10d) 


For 


^12 ® (6.8) gives A in terms of 


A = 



b 


( 6 . 11 ) 


In (I), b was assigned a value of 0.125, and therefore 
A = /2 . Table 1 presents a comparison between the results of experi- 
mental investigations of nearly homogeneous, constant shear flow tur- 
bulence by Champagne, Harris and Corrsin^^, and Harris, Grah^ and 
Corrsin^® and the nondimens ional turbulent velocity correlations pre- 
dicted by Eqs. (6.10) and Eq. (6.11) for b = 0.125. 
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Table 1 


Comparison Between Measured and Theoretical Turbulent 
Velocity Correlations for Constant Mean Shear Flow 



DATA 


THEORY 

Turbulent 

Velocity 

Correlation 

Champagne , 
et al.^’ 

Harris 
et al.^® 

(b = 0.125) 

/ 2 

u^u^/q 

0.47 

0.50 

0.50 

/ 2 

0. 25 

0.20 

0.25 

/ 2 

UjUa/q 

0.28 

0.30 

0.25 

/ 2 

UiU^/q 

-0.16 

-0.15 

-0.176 


Agreement between theory and data is acceptable for pur- present 
purposes. 

Having an approximation for the zero- separation limit we now can 
go on to select the functional representation for the two-point, 
two-time correlations. 


7. SELECTION OF TWO-POINT, SPACE TIME CORRELATION FUNCTION 

One of the basic requirements of our analysis is that we can ’ 
select a functional form for Qii ( . ^3 ; t) that is both versatile 
enough to characterize the behavior of the two-point, two-time velo- 
city correlations, yet be amenable to analysis. The trial function 
that we select should be functionally explicit in separation space 
(so that spatial integration can be accomplished analytically) and 
contain a set of separation- time-dependent parameters. The form of 
the selected function is based on examination of experimental measure- 
ments of two-point, two-time velocity correlations, supplemented by 
theoretical conditions which may be inferred from continuity, symmetry, 
homogeneity and stationarity of the turbulence and the limiting fom 
of the correlations at zero and infinite separation in space and time. 
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The method of solution then consists of forming a set of equa- 
tions based on taking the (m,k)-fqld spatial moments of the govern- 
ing equations containing the selected trial function, where m 
refers to the moments in the spatial coordinates and k the direc- 
tional dependence and then integrating in separation space to form 
equations dependent only on t . The product m x k must be chosen 
such that the nvnnber of equations formed, in conjunction with con- 
ditions derived from supplementary constraints such as continuity, 
are sufficient to calculate the total number of parameters contained 
in the trial function. 

Since in any specific turbulent flow the characteristic size of 
the energy- containing eddies is set physically, such as by the largest 
dimension of the flow, by the diameter of a jet, etc., 
we must leave one parameter of the analysis free to be specified. In 
general, this should be one of the integral scales, and this is the 
approach taken here. 

7.1. Supplementary Constraints 

The form of the two-point, two-time velocity correlation trial 
function must conform to several conditions 

1. Loss of correlation at large separation 
in space and time 


Qij (5^ ; t) ^ 0 

(7.1a) 

Qij ^ ® 

(7.1b) 


2. One-point, one-time limit 


Q^j(0 ; 0) = (u.u^)q (7.2) 

3. Incompressibility 


8Q.. 3Q.. 

^11 = _il = 0 






A. Symmetry with respect to the 
unidirectional shear flow 


(7.3a,b) 


plane for 








= 0 


(n odd) 




(7. A) 
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5. Homogeneity and stationarity of turbulence. In 
coordinates convected with the mean flow 


^ij ^ij 

6. Taylor's hypothesis; temporal history at a 

fixed point is related to the convected spatial 
structure by 



(7.6) 


The decay of Q . . with t in convected coor- 
dinates makes tMs only approximate. A memory 
function gy (t) will be used to account for 
this decay in separation time in a convected 
coordinate system. 

7.2. Previous Approximate Functions 

Frenkiel^ ® ^ ° evaluated a family of approximations for the iso- 
tropic lateral velocity correlation g(?) of the form 


g(C) = exp (7.7) 

where 

i|^(^) = a^ (m^c^) (7.8a) 

or 

= 1 +J]a^ c^Ur (7.8b) 

He found adequate agreement with wind tunnel data could be achieved 
near ? = 0 only for m > 2 , but over the whole curve an adequate 
approximation could be found for 1 < m < 2 . 

For aerodynamic sound analyses, Ribner^^,^^ approximated the 
fluctuating pressure correlation using the Gaussian form 

p2 = p2 exp [-a^(C;L - Ut)^ - a^C^ - - a^a^U^x^] (7.9) 

where provision for different turbulent scales is made in the factors 
a^ , and ^3 * 
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Several investigators have based their calculations on locally 
homogeneous isotropic turbulence correlation models. ^ ^ ^ > ^ ^ ^ ® 

Goldstein and Rosenbaum^® treated a more refined model of axisymmetric 
turbulence, using a zero time delay function similar to Eq. (7.7) with 
m = 1 and (}> (x) = 1 . 


7.3. Present Model 

The trial function selected here combines some of the features 
of the Frenkiel and Ribner models described above. The rationale 
underlying its adoption is discussed below and its limitations will 
be presented shortly. 

The model we will use here is 


Qij ^ij (7.10a) 


Rij (^1 


- a..(K. - U, - 3.. Co - Y-.C^ 

ij Vl 1„ / ij^2 


ij^3 


- - “l. 


IJ 

-2 

’2 


•exp - 2 




2 

■ • 2 ^ j ^ 


• •o’* 


ijl 


(7.10b) 


where there is no svimmation implied in Eq. (7.10a). The spread 
for the ij correlation in the kth coordinate direction is simpli- 
field somewhat by writing it 


a 


ijk 


‘^ij^k 


(7.11) 


Thus, the anisotropy in spread is assumed equal among the different 
correlations. 

The model guarantees that the correlations approach zero at 
large space separation since they decrease at least exponentially 
with • gij ("t) . the memory function, may be calculated numeri- 

cally but will be evaluated near t = 0 as we show later. Qij ap- 
proaches zero exponentially with t 2 and with gij(T) unless we are 
translating with the flow at velocity , in which case loss of 
correlation is dependent primarily on g^j (x) . The gradient of 
U(t 2 ) is neglected in the trial function for this initial study. It 
will be added later if it appears that doing so can improve the pre- 
sent model. 
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The bracketed polynomial in Eq. (7.10) permits zero crossings 
which are necessary to satisfy continuity. Equation (7.10) does not 
do this in the differential sense (Eq. (7.2a or b)) but we force con- 
tinuity in an integral sense as follows. From Eq. (7.3a) we integrate 
the equation over half- space for each coordinate direction 


00 oo 



d5 


& 



= 0 





(7.12a) 


(7.12b) 


(7.12c) 


The alternate statement of continuity, Eqs. (7.3b), is not independent 
of Eq. (7. a) (see Batchelor^ p. 27). The results of Eqs. (7.12) are 
five independent equations in the model parameters. 


The continuity constraints require that no net flew passes throu^ the 
three perpendicular coordinate planes passing through ^ = 0 . Since 
we may translate these planes freely in the homogeneous flowfield 
without effect, these conditions are met everywhere. Note that these 
conditions are x dependent. 


Qii is expressed in Eq. (7.10a) as the product of the zer o sep-; 
aration, one-point, one-time double-velocity correlation - (Uj^uj)© . 
a normalized function which we designate as Rij (^1 ^ 59 » ^3 > "r) > 

" '0,0,0;0T = 1 , and we 


the normalized memory function gij(x) . Rij (0 
specify gij(O) = 1 . 


Taylor's hypothesis is satisfied by writing the flow direction 
argument as (5l ■ ^lo^)^ • ^ii symmetric in this argument since 

T changes sign with 5i • In calculating the dependent parameters 
it is important to separate the variations that occur due to convec- ' 
tive effects, i.e., separation coordinate changes, from variations^ 
that occur in convected coordinates. It is the latter that consti- 
tute. the actual variations that are of interest. Thus, in solving 
for the time variation of the function parameters the convective 
velocity is set equal to zero. The decrease in Qj[j with convective 
separation is accounted for by the convective argument (5]_ ■ lo"^) • 

gi^(x) provides decorrelation in the coordinate system translating 
with the flow. 
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One of the Cj. functions is set equal to one both initially 
and for t > 0 . At t = 0 we specify one length scale in the same 
direction. For x > 0 the memory function is used to provide the 
remaining time dependence of the length scale in the appropriate 
direction. The spread functions are not taken to be time de- 

pendent, being set oy specification of the initial conditions. 

These initial conditions are calculated from the governing equations 
with separation time derivative equal to zero. ( Qij/9x =0 at 
T = 0 since the flow is assumed homogeneous and stationary) . Note 
that initially only one parameter is specified, the integral scale 
of turbulence for one ij correlation in one coordinate direction. 
The integral scales in the k = 1, 2 and 3 directions are defined by 
(no summation) 


00 

A^j^(x) = \j^^l “ Ux,0,0;T)d(5^ - Ux) 

00 

A^^^(x) = R^^(0,52,0;x)dC2 . 

oo 

Ry (0,0,?3;x)d?3 


(7.13a) 


(7.13b) 


(7.13c) 


(k) 

One value of Ay (0) will be specified by integrating a sel- 
ected experimental measurement of Rij(?k;0) over ^ 1 ^ . Note that 
this is not a universal specification but must be done for each 
individual flow analyzed. 


The a^. , $y and Yij parameters, which are x dependent 
serve to satisfy continuity in the integral sense discussed above. 

will have zero crossings dependent upon the determination of the 
ay , By and Yij functions, (as well as on • 

7.4. Model Limitations 


The most serious limitation of the model is its restriction to 
only the integral scales of turbulence. Thus, although the behavior 
of the energy containing eddies can be represented by the model, 
there is no explicit dependence on the turbulent dissipation scale, 
or microscale. Therefore, we must expect that the chosen double- 
velocity correlation function will not provide a good representation 
of data in the region of small spatial separation. This can be reme- 
died by adding additional terms to our model which are scaled by 

rather than by A where is the directional microscale 

T 1 •'ll 11 
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for each Qij . Indeed, the ideal representation would be an inte- 
gral function which incorporated the entire spectrum of scales which 
occur in the flow and which account for the behavior of the turbul- 
ence in wavenxmiber space. However, for this initial study the 
analysis will be based on an integral scale model only. Note that 
better agreement could be achieved near = 0 hv using an 

exp (-l^i^l) type dependence. A Gaussian was chosen here in antici- 
pation of adding the microscale dependence once our concept has been 
proven valid. 


8. MOMENTS OF GOVERNING EQUATIONS 
8.1. Convention for Nondimensionalization 

Equations (7.10) are now substituted into Eq. (5.5) and inte- 
gral moments are taken with respect to the coordinate directions. 
Before doing this, it is convenient to nondimens ionalize the vari- 
ables by a length scale q/ (dU]^/dx 2 ) , a time scale (dUj^/dx 2 )“^ , 
and a velocity scale q . 

The definition of A , Eq. (6.7) is carried forward and we 
define 

V 

N = ;^ (8.1) 

The parameters used to nondimens ionalize the variables of the 
analysis are given below in Table 2. 

With the understanding that all variables appearing henceforth 
are nondimens ionalized unless noted otherwise, the model equations 
governing Q^j appear as 


3Q.. 3Q.. T 

-^ + (U + 52> = - (6iiQ2j + «jiQi2> - AlQij - 


where 


+ NV^Q. . - 2bAQ. . 


Q.. =u..R..g..(T) 


(no sxmimation) 


and 

Rij = tl - a.j(Ci - Ut)2 - 

f E^] 

• exp I- - 2 “ i (5l - Ut) + “2 + “2 

^ij t^l 2 ^3^ 


( 8 . 2 ) 

(8.3) 


(no sxraimation) (8.4) 
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TABLE 2 


Nondimensionalization of Variables 


Variable 

Nondimens ionalized 

Nondimens ional 


by: 

notation 



Q. . 

^13 


q 

U 


2 

q 

u. . 

13 

a . . 

13 

[q/(du^, 

/dx 2 )]^ 

a . . 

13 




3. . 

ij 


N 

f 

Y • • 

13 


q/(dU^ 

/dx 2 ) 

13 





a . . 

3-3 

1 

N 

/ 


T 

(dU^/dx 2 )'-^ 

X 


40 




8.2. Sxjmmation Convention 

Since it will be convenient from this point on to discard the 
summation convention when denoting the velocity correlation sub- 
scripts i and j and the direction subscript k , these will be 
reserved exclusively for these purposes. Therefore^ whenever i , 
j and k appear as subscripts, no summation will be implied. 

Other subscripts denote the svimmation convention as usual. 

8.3. Moment Function 

There are no set rules regarding which moments will give an 
optimum result when using this method to obtain solutions to any 
given problem. Given the number of unknowns and the governing dif- 
ferential equations determines the number of moments that must be 
taken to obtain a sufficient set for solution. The philosophy used 
during this analysis was to develop a family of rational moment 
functions and to avoid a partial application of any subset of this 
function. For instance, if the function were ^ ^ ^ would have 

to range from 1 to 3, and not just from 1 to 2, etc. 

The function chosen can be written as n(mj^) where the argument 
denotes the product of coordinates of subscript k each to the power 
m^ , k ranging from 1 to 3. 

n(mj^) = n(m^,m2,m2) 

thus for example 

IT(0,1,2) 


The integrated moment of Q^j with weighting function n(mj^) is 
then defined as 

00 OO CX) 

h? “/// "‘V'5i3'i5l‘l52d53 (8.6) 

— 00 — 00—00 

2 

The integral moment of Q^j with weighting function, ^2^3 
therefore expressed '• 

^ij 

Other integral moments that will be required are defined as 

OO OO OO 

(tn, ) f r f 9Q. 

Ci/ =Jjj I'V'" + Tjf '8'” 




(8.5) 


• 0.1.2 


^2^ 


3 = ^2^5 
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00 CO oo 





( 8 . 8 ) 


8.5. Family of Moment Equations 

Taking the moments of Eqs. (8.2) with respect to the weighting 
function II (mj^) then leads to a family of equations of the form 



(nik) (mi,^ 

il 2j jl i2 


r (ir 


<’"k> U 

3 ij US, 


+ 


ND 


<"k> 

ij 


(tUk) 

2bAI . . 


(8.9) 


(nik^ 

Integration by parts may be used to evaluate , while 

Green's theorem is used to determine Dj^j rather simply. Over 
volume V 


///[ n(mj^)v2Q^j - Q.jV^ n(mj^)jd?^dC2d?3 

V 

= n ” ^ij'^ H(nij^)|dS 


( 8 . 10 ) 


where S is the area bounding the volume of integration V and n 
is the outward pointing unit normal to S . Letting V go to in- 
finity, the bounding integrals vanish since Q^. approaches zero 
exponentially. Then 


00 00 00 


/// n(mk)V^Qijd5idE2d53 

OO 00 00 


— oo — 00—00 oo oo oo 


( 8 . 11 ) 


-00 — oo -00 
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Starting with “ (0>0»0) the moments can be evaluated in 

iiicr easing order. 


The lowest or zero moment mk = (0,0,0) is related to the 
anisotropic scales > where the moments are taken along the 

coordinates, see Eqs. '^(7.12). The first moments in each coordinate 
direction mk = (1,0,0) 5 (0,1,0) ; (0,0,1) are zero since the tur- 
bulence is homogeneous and these are odd moments. The first mixed 
moments, mk = (1,1,0) ; (1,0,1) ; (0,1,1) make a contribution due 
to the skew symmetry produced by shear in the (Cl»C2) plane (mani- 
fested in the correlation function by the presence of V>ij). The 
family of moments required for this analysis ends with the second 
moments mk = (2,0,0) ; (0,2,0) 5 (0,0,2) giving a sufficient set 
of equations to provide solutions for all parameters. 

Before evaluating the integrals defining the moment functions 
some definitions are made for convenience. Let 


d . . 

ij 

^2 2 
ij 1 

(8.12a) 

Ak 

e. . 

2 2 
ij a-J 2 

(8.12b) 

Y • • 

3-J 

2 2 
a-j ij 3 

(8.12c) 

"ij 

^2 

IJi • • CT • • Cl Co 

a.j 1 2 

(8.12d) 


The correlation function, Eq. (8.3), is now substituted into : 

the integrals defining and to relate 

them to the correlation function parameters. Note that these inte- 
grals are to be evaluated for U = 0 since the parameters are being 
calculated in a convected coordinate system moving with the mean 
flow. 


mi, = (0,0,0) - 

anisotropic scales 


Th 

A 


l^pecial 

ij 


ized form of the zero moments are the 



(8.13b) 


A<?> 

ij 




A<?> 

3-J 


r o -1 

■ 5 - C 3 °ijI 2 - Yyl 


( 8 , 


~ (1>0»0) ; (0,1,0) ; (0,0,1) - As noted above the first 
moments are zero. 

^^k = (1,1,0) ; (1,0>1) ; ( 0 , 1 , 1 ) - Only the first of these 
three mixed moments is non-zero and results in 


( 1 , 1 , 0 ) ^3/2 3 

^ij “5 ^ij^ij^ij^l^2^3^ij 


(8. 


(nik) (m, ) 

The associated values of C.. and D.. are 

iJ ij 

(1’1’0)_ 7^3/2 3 3 

^ij 2 ^ij^ij^ij‘^1^2^3 ■ 2^“ij ^^ij 


(1,1,0) 

D. . = 0 

ij 


( 8 . 


my = (2,0,0) ; (0,2,0) ; (0,0,2) 


(2,0,0) _ 7j3/2 

^ij " ~T~ 


o ^ 

0^0203 [1 - + Y^*)] ( 8 , 


o V P . . Y . • 

2" ij ij ’rj 


(0,2,0) 3/2 3 

^ij " ^ ^ij^ij°ij''l"2''3 


^ ^ii^ii°ii^l^2‘^3 " 2^“ij ■*■ ^^ij 


(0,0,2) 3/2 




7T 


^ij®ij"^ij ^ 1 ^ 2^3 ~ 8 ,. ^ + 3y,. ^)] (8 


2v-ij - 




13c) 

14) 

15) 

16) 

17a) 

17b) 

17c) 
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(8.18a) 


( 2 , 0 , 0 ) ^ 3/2 5 


( 0 , 2 , 0 ) ( 0 , 0 , 2 ) 


(8.18b,c) 


( 2 , 0 , 0 ) ( 0 , 2 , 0 ) ( 0 , 0 , 2 ) 

D. . = D. . = D. . 

ij ^3 iJ 

Q/9 Q ^vi 


(8 . 19a,b , c) 


The mixed moment (1,1,0) provides one set of equations in the 
11,22,33,12 components. The second moments provide three equations 
for each ij component corresponding to the three coordinate direc- 
tional moments. Note that q^ 2,0,0) _ and 

_ j( 0,2,0) ^ coupling the mixed and second moments. 

8.6. Integral Continuity Constraints 

Substitution of the correlation function, Eq. (8.3) and Eq. (8.4), 
into the integral constraints, Eq. (7.12), yields the following five 
equations 


hi 

+ Yii = 2 

(8.20a) 

^22 

+ Y 22 = 2 

(8.20b) 

^^33 

+ ^33 " ^ 

(8.20c) 

“^12 

+ Yi2 = 2 

(8.20d) 

&12 

+ Yi2 = 2 

(8.20e) 


8.7. Method of Solution 

The integral moment equations , the continuity constraints and . 
only one of the integral scales are sufficient to define the correla 
tion function parameters by integrating Eqs. (8.9). At t = 0 , 

3 Q. /3x = 0 is the initial condition required. Differentiation of 
the^ correlation function with respect to t yields the conditions 
that a..(o) =0 , B.-(o) =0 , etc., for this initial condition to 
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be valid universally. With these conditions Eqs. (8.9) become a set 
of nonlinear algebraic equations in the correlation parameters at 
T = 0 which may be solved to determine their initial values. If it 
is assumed that the acoustic sources are compact then a solution 
derived by expansion of the variables about t = 0 can be determined, 
which has the advantage of providing an analytical solution for the 
higher derivatives in t . Since the fourth derivative with respect 
to T is required to evaluate the acoustic integral, this approach 
will avoid nxraierical evaluation of the derivatives. 


Since the first derivatives are zero at x = 0 the correlation 
parameters which are chosen to be functions of time in this simpli- 
fied analysis are expanded in series in x as follows 







a . . 



1 


+ 

+ 

+ 

+ 


/s 



2 

2 

2 

2 


+ 

+ 

+ 

+ 





4 

4 

4 

4 


+ * • • 

+ . . . 

+ . . . 

+ . . . 


(8.21a) 
(8.21b) 
(8. 21c) 
(8.21d) 


The decorrelation function b is x dependent, as it must be 
for giT to vary from its initial value. Expressing b as a 
series in x/xq , where xp is the characteristic turbulence time 
A/q , the nondimens ional expansion is 

b = b^ + b^Ax + b2A^x^ + b^A^x^ + b^A^x^ + ... (8.22) 

where b^ are specified constants. Nothing is known about this 
function at the present time and definition of the physics and form 
of this parameter is put off to a following study, bj^ will be 
chosen by matching the results of the present analysis to the distri- 
bution of measured acoustic data. 

Substitution of the expansions, Eq. (8.21) and Eq. (8.22), into 
the moment functions, Eq. (8.14) to Eq. (8.19), yields series 

expressions for and in powers of x with 

coefficients which are functions of , etc. 
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(8.23a) 


/“k> 


(nik) (V 2 , j^^k\3 , . 

I.. +1.. T +1.. T +1.. T + 


3-J, 


1J2 


XJ3 




IJ 


(tn. ) 2 3 L 

C..^ + C.. + C..*^ T-^ + C..^ + 


3-J, 


ij 


IJ3 




(8.23b) 


(tn, ) 
D..*^ 
iJ 


(nik) (m, ) 2 (mi.) 3 (m, ) , 


^3, 


3-3 2 


133 


^Ja 


(8. 23c) 


Separation of Eqs. (8.9) into terms in powers of t then pro- 
vides a set which may be solved for successively higher coefficients 


o 

T : 


(nik) 

C. + 
3-3^ 


, <V . , <“k>i 


+ A 


T^“k> 1, /“k^- 


(mk) (m, ) 

-ND..‘^-K2bAI..*^ =0 

3-3^ o rj^ 


(8. 2Aa) 


<”k> 

3-3 2 


9 (mi-) 
b,A^I. 

1 3.1 


(8.2Ab) 


(nik) 

31. 

3-J3 


= - C 


3-3 2 


(mk) 

6 . I Q . + 

il 2 j2 


(m, ) • r (m. ) , 

-A I. .‘^-46.. I 

3 1 3.22 13 2 3 ij 


(m^)l 


(nik) 

-1- ND . . - 2A 

3-3 2 


(m, ) 2 ^™k^l 

b_I..^ +b,A^I..‘^ 

^3o . 


i/ ^ • 

o 1 J 2 


(8.2Ac) 
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1J3 


<”k> 


<”k) 


+ «ji^i2; ] - A 


f <”k> u 'V' 
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+ m.. -2k 

XJ3 
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3 (mi,)' 
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o 3.J3 1 3.J2 3 

(8.2Ad) 
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Solving the nonlinear algebraic set of zero-order equations, 

Eqs. (8.24a), provides the initial conditions for the correlation 
parameters at t = 0 . The model constants that must be specified 
are A , bo , u^^ and v^ . The solution provides ai-i , Bi-] , 

^ijo » ^ijo ’ ^ijo ^^o * initial value of is 

specified as gij (o) = Sijo “ For now, we assume 

and ck do not vary with Separation time, affording considerable 

simplification without undue loss in generality. 

Since the higher-order coefficient equations are linear in those 
of lower order, Eqs. (8.24) provide the series expansion coefficients 
needed to define the separation- time derivatives at t = 0 using 
elementary matrix algebra. 

9. RESULTS AND COMPARISON WITH TURBULENT SHEAR FLOW DATA 

Tables 3 and 4 present the zero- order coefficient solutions to 
Eq. (8.24a) for model constant values A = /2 , un = 1/2 U22 = 1/^ i 

^33 = 1/4 , u12=-1/(4/2) , vc = 0.3 and three values of bo: 0.0, 
0.02 and 0.05. Figures 3 and 4 are plots of the correlation function 
Rij for bo = 0.0 and 0.05, respectively. Comparison of these 
figures with measured two-point, one-time turbulence correlations 
makes possible an evaluation of the qualitative aspects of the re- 
sults. Such measurements in a two-dimensional constant shear flow 
have been made and reported by Champagne, Harris and Corrsin^^, 

Harris, Graham and Corrsin^®, and Rose^®. Figure 5 presents mea- 
sured values of Rn and Rl 2 measured along the coordinate direc- 
tions at T = 0 as measured by Harris, et al.^® The directional 
behavior of these correlations are typical of those measured by 
Champagne, et al.^^ and Rose et al.^ The measurements in the ^3 
directions, normal to the mean flow and velocity gradients, decrease 
most rapidly with separation distance and exhibit extensive regions 
of negative correlation. Those in the flow direction, , are 
sloxvest to decorrelate. These characteristics are present in gen- 
eral in the theoretical predictions, although the trends are some- 
what exaggerated. Since Rij(t;0) < Rij(0,0) in a physical flow, 
the theoretical behavior of Rs 3 ( 0 ,0 ,C 3 ; 0) in both Figures 3 and 4 
is the most serious discrepancy in the model. It should be remem- 
bered, however, that the present model is an integral representation 
and as such may exhibit local inconsistencies with little serious 
effect in its portrayal of overall behavior. Including the effects 
of a turbulent microscale explicitly in the correlation function 
would probably eliminate the undesirable overshoot in R 33 ( 0 , 0 , C 3 ; 0) . 
Examination of the model also reveals that employing a directionally 
dependent value of v^ could reduce the negative overshoot in 
F-ii ( 0 , 0 , ? 3 ; 0 ) while increasing the rate of decorrelation in 
Rij ( 51 , 0 , 0 ; 0 ) with separation distance, thereby improving comparison 
between measurement and theory. This was not done for the present 
until more empirical evidence can be found to justify this extension. 
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Figure 3. Theoretical two-point correlations at t = 0 for ij = 11,22,33,12 and 
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Figure 5. Measured two-point correlations at 
from Harris , Graham and Corrsin. 
q = 90.5 cm/sec. dU,/dx 2 = 44 sec 




TABLE 3 


Trial-Function Coefficients at Zero Time Separation 


ij 

a . . 

^Jo 

0.. 

^Jo 

/V 

Y • • 

’ ^ 



a • . 




b = 0.0 



11 

0.1227 

0.1169 

1.8831 

-0.0990 

1.1985 

22 

0.1746 

0.1151 

1.8254 

-0.1495 

1.1675 

33 

1.308 

0.6916 

-1.9082 

-0.2216 

1.0254 

12 

0.1186 

0.1186 

1.8814 

-0.1028 

1.1959 




b = 0.02 



11 

0.1068 

0.1028 

1.8972 

-0.0784 

1.1896 

22 

0.1500 

0.1015 

1.8500 

-0.1086 

1.1687 

33 

1.3277 

0.6723 

-1.7699 

-0.1730 

0.9739 

12 

0.0037 

0.0737 

1.6768 

-0.0237 

1.:|.882 




b = 0.05 



11 

0.0743 

0.0721 

1.9279 

-0.0386 

1.1706 

22 

0.1006 

0.0718 

1.8994 

-0.0493 

1.1624 

33 

1.3726 

0.6274 

-1.5060 

-0.0572 

0.8522 

12 

0.0728 

0.0728 

1.9272 

-0.0382 

1.1700 


& 
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TABLE 4 


Scale Factors at Zero Time Separation 


b(o) 

k 

0.0 

0.02 

0.05 

1 

3.6313 

4.4728 

6.8162 

2 

1.0000 

1.0000 

1.0000 

3 

0.3028 

0.2830 

0.2358 


10. RADIATED ACOUSTIC POWER BASED ON 
TURBULENT CORRELATION FUNCTION 

From Eq. (3.8)' the acoustic power radiated in direction (ijj.tl)) 
from a unit volxime element at y is , in the Proudman formulation 
and in dimensional form 




P, 


16 tt^c^x^ 

o 



2 ,2 
V V 
X X 


dt 


( 10 . 1 ) 


where Vx and vi are the components of total velocity at y' ^ 
and y” in the direction of ^ , see Figure 1. Since S and r 
are nearly parallel for x in the acoustic far fiel^, we shall^ 
approximjite X by r . Let the total velocity at 5 = 0 be v 
and at 5 , . The mean flow is along 5l so that the compon- 

ents of total velocity are 


u. + 

(10.2a) 

“i + «lk«' 

(10.2b) 
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where and are the fluctuating components. The velocity 
components and Vx in terms of vj[ and are, to a first 
approximation*'' 


V = V-, cos \p 
X 1 ^ 

+ V2 

sin 


cos cj) 

T v^ sin 


sin 

«}> 

(10.3a) 

v' = V, cos 

X d. 

+ V' 

sin 


cos (|) 

+ v^ sin 


sin 

<!> 

(10.3b) 


2 2 

The correlation v^Vx is found by substituting Eq. (10.2) into 
Eq. (10.3); squaring Eq. (10.3a) and Eq. (10.3b) and forming their 
product and taking the e nsemble av erage . Th e result is ano express 
sion containing terms in u,.u^u!uj|, , UU'u.uI , U^U'^ , U^u. , U2ul 


UU^Uj 


71 


and 


1 P J q 1 J 

The first two of these will contribute to 


the flow noise. The following three are constant with x and will 
vanish upon differentiation. The last two terms are triple-velocity 
correlations. For homogeneous, stationary, turbulence, these are 
odd functions by invariance with reflection of coordinates. Change 
of sign of coordinates must be accompanied by change of sign of 
separation time. Since we have assxmed the correlation function 
may be expressed as the product of separation space and separation 
time variables , the triple-velocity correlations will integrate to 
zero over sepui^ution space. 

In order to provide a tractable analysis for the fourth-order 
correlations, normal joint probability of ui and uj is assumed. 
This is done based on Batchelors^’ argument that the part of the 
joint probability distribution of the velocities associated with 
the energy containing eddies is approximately normal at a fixed time 
and at points sufficiently separated in space. This assumption has 
been used by numerous invest igators . ^ ® ^ ^ ^ ® ^ ® ^ ^ ® Goldstein 
and Rosenbaim^® extended this argument to time separation by arguing 
that the correlations will be subject to even more random influences 
from the neighboring flow when they are separated in time. Thus, 
their joint probability should be even closer to normal. 

Using the assumption of normal joint probability the fourth- 
order velocity correlations may be written (see Batchelor^’) 


The effect of % on the components of v. and vl in the direction 
of X is of order C/x << 1 and will bi neglected in this analysis. 
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( 10 . 4 ) 




= Q. (0)0. (0) + Q. .Q + Q. Q . 
^3-p ^jq ij pq ^iq pj 


The first term on the right will vanish upon differentiation with 
respect to t . 

T 

The other nonvanishing term is of the form UU'uiUj , where 
U' = U + ^2 nondimensionally. For simplicity the variation in U 
with is dropped and only the leading term in is retained. 

This leaves terms of the form 

i j 


The noise associated with the fourth-order turbulent correla- 
tions has been desi gnated self-noise, while that generated by terms 

of the form U^uiUj is known as shear noise. These contributions 
will be written separately in what follows. 

P('l^f<t>.y) > the sound power radiated along x in the direction 
(4),c()) per unit volxme of shear flow per unit area at x may be 
expressed as 


, L dU, 
P(1.,*.y) - ^^^2 dX 2 


ill! 


“m J 


u2(5)5^AijUi. ^ (g..Wy> 

ij 


^^ijpq^ij^pq 77^ (^ijSq^ijpq)l 

ijpq 


(10.5) 


Although not explicitly stated, u^j and Upq are functions 
of absolute space when a real flow is treated using the ''locally 
homogeneous" approximation. 


Aj[j an d B^^ pg are functions of ip and <() as derived by ex- 
pansion of v^v'^ The sxmmation sign indicates summation over all 
correlations that contribute to the sounu power, P(if',<}>,y) is non- 
dimens ionalized as follows 


P(4',4>.y) 


(»,♦,?) 

P u® 
o m 


( 10 . 6 ) 
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and ^ijpq integrals of and over sep- 

aration space 


W. . 

ijpq 


00 00 00 

'/// 

— 00 — 00—00 
00 00 oo 





(10.7a) 

(10.7b) 


Ujh is a reference velocity and Ljjj a reference length. Neither of 
these is a function of y . 

9 9 

The values of Aij and derived by expanding VxVv 

the correlations which have been"' included in this analysis and which 
contribute to the sound power are shown in Tables 5a and 5b. The co- 
efficients in these tables include the permutations of the indices of 
the correlations to include all contributing terms. 

If the shear flow is axisymmetric the sound power may be aver- 
aged over <|) , since the nonaxisymmetric power from individual vol^Jme 
elements of the flow will mutually cancel on a time average basis. 
Define 


2TT 


A. . = 


2tt J 


2tt 


®ijpq - L / ®ij 


pq 


d(j) 


(10.8a) 


(10.8b) 


The results of these averages are shown in Tables 6a and 6b. The 
(})- averaged sound power is expressed as 


^'“ 117^34 ^ (g,jW 

Z \jpq“ij“pq ^ (®ij®pq”ijpq)} 

-1 T'rvn O ^ ^ 


P('l^, 




(10.9) 


ijpq 
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Sound Directivity of ij and ijpq Components 


TABLE 5a 


ij 




A. . 

3-J 


11 

4 

4 

cos 




22 

4 

cos^ 


2 2 
sin CO s 

4> 

33 

4 

cos^ 


sin tj; sin^ 


12 

8 

cos^ 


sin cos (f) 



TABLE 5b 


ijP.q ijpq 


1111 

2 

cos^ 






2222 

2 

sin^ 

'P 

cos^ 

4> 



3333 

2 

sin^ 

'P 

sin^ 

■<P 



1212 

8 

cos^ 


sin^ 

<P 

cos^ 

<P 

1122 

4 

cos^ 


sin^ 

4> 

cos^ 

<P 

1133 

4 

sin^ 


2 

cos 


sin^ 

4> 

2233 

4 

sin^ 


sin^ 

<P 

2 

cos 

4> 

1112 

8 

sin 1 

P > 

cos^ 


cos <j) 


2212 

8 

sin^ 


cos 

ip 

cos^ 

4> 

3312 

8 

sin^ 

'P 

cos 


sin^ 

4> cos 
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(})-averaged Sound Directivity of ij and ijpq Components 


TABLE 6a 


ij 



11 

4 

4 

cos 




22 

2 

2 

cos 

'p 

sin^ 


33 

2 

2 

cos 


. 2 
sin 

'P 

12 

0 






TABLE 6b 


ijpq 


R. . 

ijpq 


1111 

2 

4 

cos 



2222 

3/4 sin^ 


3333 

3/4 sin^ 


1212 

4 

cos^ 


sin 

1122 

2 

cos^ 

’I' 

sin 

1133 

2 

cos^ 


sin 

2233 

% 

sin^ 



1112 

0 




2212 

0 




3312 

0 
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Integration of R. . and R-^R„a separation volime 

yields ^ 


^ij " " 2^°^ij ■*■ ^ij ■'■ ^ij^^ 


(10.10a) 
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ijpq 
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/s /\ 
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(10.10b) 


where 


2 

K . . 

ijpq 


^2 „2 

a . . 0 

IT pq 


2 . „2 
o . . + o 

, ij pq 


(10.11) 


The fourth- derivative of these functions at x — 0 are obtained 
by substitution of the series expansions for the correlation para- 
meters, Eqs. (8.21), taking the derivatives and setting x = 0 . 

The results of these operations are 


^ (g.jW..) = 24H3/2c,C2C3O^gy^s(0>S<«)j 




(10.12a) 
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,3/2. 


W.. \ = 24it'" c«Co<’. . _ 

xjpq^ 12 3 ijpq 
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1-2 
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2 pq \ 


G> + <3-G> + <Y^j 
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A A 


A A 
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+ <a y--G> + <3.-Y„„G> + <Y,. i3„_G> 
pq XJ XJ ' pq XJ pq 


(10.12b) 

where the various functions in these expressions are defined by 
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^34 
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pq 
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■ [* ■ K”«o 

■-K' 

■1 
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1 


-1 


-2 
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pq' xj' " XJ pq 


(10.13a) 

(10.13b) 

(10.13c) 

(10.13d) 

(10.13e) 


<Xi-G> = X,.^G^ + X,.^G2 + X,.^ 


(10.13f) 


<Y G> = Y G/ + Y„„ G« + Y„^ 

pq pqo ^ pq2 ^ P'14 


(10.13g) 
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<X.,Y„^G> = X.. 

13 pq ijo 


Y G, +( X.. Y„^ + X.. Y^^ \G, 

Pqo ^ V ^3o Pq2j ^ 


+ X. . Y + X. . Y + X. . Y„„ 

IJ 2 Pq2 P^4 


(10.13h) 
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G. = g., . + g.. g. 


+ Sr 


4 °iJ4 ■ °ij2°Pq2 ' °P^4 
takes on the values of the expansion coefficients a 


X. . 
ij .. 

and" Y^j where appropriate, while Yp^ symbolizes 


ij 


(10.13j) 

, 3, 
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^n 


pq 
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13 
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pq 


n 


and Y 


pq 
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The sound power intensity I(ip) for axisynraietric flow may now 
be found by integrating Eq. Q0.9) over the voltme of absolute space 
containing the turbulence. Since u^j, , Upq and may be func- 

tions of f for "locally homogeneous’' turbulence, it is more con- 
venient to nondimensionalize the one-point, one-time correlations by 
the reference mean velocity Ujji . Thus, can be written 




1 Po^m 

1 0 m 

fff 


dU(y)\ 


16.2 cVl^ jjj 

''f 
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q(y)V 
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ij 
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^1 

Mm ) 2 

^ ijpq [ 

'^m 

j 


r. .W. .) 

*13 ^2 


(10.14) 


^SijSpq^ijpq)/ dy 


i3pq 


where the integral is taken over the flow volume Vf and (u^uj)^ , 
(^(yF)^ , q and the gradient dU^/dx 2 may be functions of y . 

This completes the equations necessary to compute the sound 
power, given the one-point, one-time turbulence velocity correlations 
and the mean shear flow distribution. In the next section the theory 
is applied to a shear flow and the acoustic power emission is 
calculated. 
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11. PREDICTED ACOUSTIC POWER FOR SIMPLE SHEAR FLOW 

The sound power intensity equation, Eq. (10.14), requires the 
spatial distribution of turbulence and shear. The test case selected 
to evaluate the model developed here is an axisymmetric , annular shear 
layer. The three-dimensional turbulence in this annular layer was 
calculated based on a one-dimensional shear. Figures 6 and 7 illu- 
strate the geometrical details. Figure 6 shows a one- dimensional 
velocity profile schematically with boundary conditions U(-<») = 0 
and U(+») = U . The turbulence was calculated using the second-order 
closure model of Donaldson, with the mean profile allowed to evolve 
to a self-similar distribution. The shear layer was treated as axi- 
symmetric as shown in Figure 7, where the annular layer is treated as 
thin, Ag << R , thus justifying the one -dimensional mean profile 
calculation. The core velocity is uniform and equal to Uc • 

The resulting turbulence and mean velocity gradient was then used 
to calculate ICil') as given by Eq. (10.14), with Um/Uc = % . Since 
the T dependence of the decorrelation function b is unknown, the 
test case provided the opportunity to empirically assign values to 
the expansion coefficients b^ . This was done by choosing those 
values of the coefficients which provided agreement between the cal- 
culations and an actual sound power intensity measurement in decibels. 
The test data selected were those reported by Lush. These measure- 
ments were carried out using an axisymmetric subsonic jet as the 
sound source. The annular shear layer assimiption of the present 
analysis precluded a direct comparison of the actual magnitude of 
sound intensity. Therefore, the coefficients were assigned to pro- 
vide agreement with the measured 4^ dependent directivity. This was 
done by referencing the intensity to a selected value at = 0 

dB(i|») - dB(o) = 10 (11.1) 

Figure 8 presents the comparison obtained for the following b(x) 
expansion near t = 0 

b(T) = 0.02 [1 + 16At - 6A^t^ - 2A^t^ ] (11.2) 

i.e., bo = 0.02, and specifying dB(o) = 100 for Vj , the measured 
jet velocity, equal to 300 m/sec. The predicted distributions for 
V • = 195 and 125 m/sec. were obtained by decreasing the intensity using the 
ritio of jet velocities to the eighth power. Note that the predic- 
tions contain no convection or refraction effects. 

The separate contributions of the self noise and shear noise 
may be seen in Figure 9. For the values of b^ selected to match 
the measured directivity the self noise is nine times as large as the 
shear noise along the jet axis. Although the total sound power dis- 
tribution is in agreement with measured values, this ratio and the 
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Figure 7. Configuration of axisymmetric flow used for 

acoustic calculations, with the assumption 

that A << Rt . Flow is laminar in the core 
s J 

region and outside the jet 
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relative contributions of the self and shear noise predicted here 
are only a preliminary assessment. Several factors which have a 
large influence on the results require further study. The Q 13 and 
Q 23 correlations have not been included at this point and the value 
of b = 0.02 was chosen arbitrarily. Further effort is required to 
define the behavior of b(T) . It should be noted, however, that 
there appears to be a unique set of coefficients in the expansion of 
b(T) which guarantees a decay in gij ( t) (i.e. , the second derivative 
of gii(T) less tha^zero) for a given choice for b(T) and the 
selected value of Aiy > which also correctly predicts the measured 
directivity distribution. The nondimens ional length scale 

a}?^( 0) = 1.0 , calculated using the data of Champagne, et al.^^ was 
u4ed. 

Additionally, the spatial integrals over the jet volime are 
based on an approximate velocity profile. Finally, as shown by 
Goldstein and Rosenbaum,^® the ratios of the correlation scales have 
a profound effect on the distribution of sound power. The values of 
the scales for the present results appear to be too large in the flow 
direction and too small in the direction perpendicular to the mean 
velocity and the velocity gradient. No information on the integral 
scales exists for the Lush tests and, therefore, the use of the non- 

dimensional scale length =1.0 was continued. Investigation 

of the decorrelation function b(x) and the influence of the ani- 
sotropic scales in the generation of the sound power should be part 
of the next phase of study. 

12. CONCLUSIONS AND RECOMMENDATIONS 

The feasibility of computing aerodynamic sound using a new 
approach for the prediction of Qij , the two -point , two-time velo- 
city correlations, has been demonstrated for a certain class of flows. 
The agreement between measured and theoretical sound ppwer emission^ 
directivity and the ability of the technique to predict the behaviot 
of Qij in spatial separation indicates that one-point, one-time tur- 
bulence models can be successfully extended to the two-point, two- ^ 
time problem. The chosen form of the correlation function is vindi- 
cated by the favorable comparison between theoretical results and 
measured turbulent correlations in three directions in separation 
space. The results provide confidence that the present approach is 
correct, although further effort is required. Three areas which are 
recommended as subjects for the next phase of study are; extension 
of the theory to noncompact sound generation; investigation of the 
decorrelation function to provide a_ physical basis for the specified 
variation of its separation- time behavior; and further consideration 
of anisotropic scale length effects. An axis ymme trie counterpart to 
the one-dimensional shear layer calculation presented in this report 
should be developed to permit application of the technique to axi- 
symmetric jets. 
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